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Effects  of  Seasonal  Land  Surface  Conditions  on  Hydrometeorological 
Dynamics  in  Southwestern  North  America 

ARO  Tracking  Number  53305-EV-PCS 

Enrique  R.  Vivoni,  Nicole  A.  Pierini,  Adam  Schreiner-McGraw,  Cody  A.  Anderson, 

Ryan  C.  Templeton  and  Luis  A.  Mendez-Barroso 
School  of  Earth  and  Space  Exploration 
School  of  Sustainable  Engineering  and  the  Built  Environment 
Arizona  State  University 

The  Earth’s  land  surface  experiences  seasonal  changes  that  affect  its  physical,  chemical 
and  biological  properties.  In  arid  and  semiarid  regions,  seasonal  variations  of  water  availability 
can  dramatically  impact  land  surface  conditions  (vegetation,  soil  moisture,  surface  temperature, 
heat  fluxes)  that  potentially  alter  vehicle  mobility,  airborne  surveillance,  target  acquisition  and 
aircraft  operations.  Our  study  focused  on  gathering  observational  data  sets  in  two  major  deserts 
of  North  America  that  will  pennit  quantifying  the  seasonal  transitions  in  land  surface 
characteristics  and  enable  the  development  of  model  applications  that  take  these  seasonal  effects 
into  consideration  for  improved  predictions.  We  focused  on  arid  and  semiarid  areas  in 
southwestern  North  America  since  these  are  water-limited  and  experience  extreme  events 
(flooding,  drought)  and  a  wet  summer  monsoon  season  from  July  to  September. 

Our  project  combined  three  approaches  for  characterizing  and  predicting  local  to  regional 
scale  impacts  of  seasonal  changes  in  land  surface  conditions.  The  first  approach  was  through  the 
establishment  of  a  dense  environmental  sensor  network  at  two  small  watersheds  in  the  major 
deserts  -  a  Sonoran  Desert  savanna  and  a  Chihuahuan  Desert  shrubland.  Each  sensor  network 
was  designed  to  accomplish  various  objectives:  1)  measure  land-atmosphere  fluxes  at  an  eddy 
covariance  site  and  its  surrounding  footprint,  2)  characterize  the  spatial  distribution  of  soil 
temperature  and  moisture  in  the  watershed,  3)  measure  the  channel  runoff  response  at  internal 
locations  and  the  outlet  and  4)  estimate  soil  moisture  at  various  levels  of  aggregation  from  point 
scales  to  the  larger  footprint.  For  each  case,  we  selected  sites  in  experimental  rangelands  where 
there  have  been  prior  observations  and  where  a  transition  in  vegetation  has  been  documented, 
typically  from  desert  grasslands  to  woody  plant  cover  (shrubs  or  trees).  The  experimental 
watershed  have  been  monitored  continuously  at  a  high  spatiotemporal  resolution  for  several 
years  (5  years  at  the  Chihuahuan  desert  scrubland  and  4  years  at  the  Sonoran  desert  savanna). 

The  second  approach  was  through  the  use  of  remote  sensing  observations  from  piloted 
and  unmanned  aerial  vehicle  (UAV)  platfonns.  At  each  site,  orthophotographs  and/or  Light 
Detection  and  Ranging  (LiDAR)  were  obtained  from  manned  aircraft  or  UAVs  and  served  to 
provide  information  on  the  site  topography,  vegetation  cover  and  species  distribution,  and 
derived  products  such  as  stream  hydrography.  Fig.  1  illustrates  an  example  from  the  Jornada 
Experimental  Range,  where  the  small  watershed  boundary  was  derived  from  a  digital  elevation 
model  obtained  through  processing  UAV-based  imagery.  Also  depicted  are  the  locations  and 
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characteristics  of  the  environmental  sensor  network  instruments.  The  local  scale  observations 
from  aircraft  or  UAVs  were  complemented  with  remote  sensing  from  satellite  platforms  that 
provide  regional  estimates  of  surface  conditions  such  as  vegetation  greenness,  land  temperature 
and  surface  albedo  that  are  useful  for  modeling  purposes. 

The  third  approach  was  through  the  use  of  a  high  performance  computing  (HPC)-based 
hydrologic  model,  the  TIN-based  Real-time  Integrated  Basin  Simulator  (tRIBS).  This  physically- 
based,  distributed  model  utilizes  triangulated  irregular  networks  (TINs)  to  describe  a  watershed 
domain  and  can  be  applied  to  remote  areas  in  arid  and  semiarid  regions  for  providing  hydrologic 
predictions.  To  gain  confidence  in  the  model,  we  tested  the  predictions  with  the  dense 
observations  obtained  from  the  environmental  sensor  network.  Fig.  2  illustrates  the  type  of 
datasets  that  are  useful  for  hydrological  model  testing  from  the  Santa  Rita  Experimental  Range. 
Model  simulations  have  been  conducted  at  two  small  watersheds,  where  the  simulated  soil 
moisture  and  latent  heat  flux  matched  well  with  the  site  data.  The  model  results  were  used  to 
illustrate  differences  in  runoff  response  due  to  woody-plant  encroachment.  Additional  model 
testing  occurred  at  the  Jornada  Experimental  Range.  Watershed-scale  simulations  that  are 
underway  utilize  the  individual  soil  profile  and  channel  sampling  locations,  as  well  as  the  basin- 
averaged  observations,  for  model  testing  purposes.  In  additions,  innovations  in  using  the  high- 
performance  model  have  been  made  for  larger  domains  where  the  seasonal  variability  in  land 
surface  conditions  was  taken  fully  into  account.  Fig.  3  illustrates  comparison  of  simulated  land 
surface  fields  from  the  physically-based  model  in  a  large  watershed  (-3500  km2)  as  compared  to 
remotely-sensed  (independent)  observations  of  these  conditions. 

Our  efforts  in  the  project  resulted  in  the  following  peer-reviewed  publications  to  major 
international  journals.  These  publications  involved  the  MS  students  (Ryan  Templeton,  Nicole 
Pierini,  Cody  Anderson)  and  PhD  students  (Adam  Schreiner-McGraw,  Luis  Mendez-Barroso, 
Nicole  Pierini)  who  were  trained  through  the  project,  along  with  collaborators  identified  as  the 
projects  evolved.  Two  additional  publications  are  currently  in  review  as  well,  as  noted  below,  for 
a  total  of  1 1  peer-reviewed  products  from  the  overall  effort,  to  date. 

Vivoni,  E.R.,  Rodriguez,  J.C.  and  Watts,  C.J.  2010.  On  the  Spatiotemporal  Variability  of 
Soil  Moisture  and  Evapotranspiration  in  a  Mountainous  Basin  within  the  North  American 
Monsoon  Region.  Water  Resources  Research.  46:  W02509,  doi:  10.1029/2009WR008240. 

Vivoni,  E.R.  2012.  Diagnosing  Seasonal  Vegetation  Impacts  on  Evapotranspiration  and 
its  Partitioning  at  the  Catchment  Scale  during  SMEX04-NAME.  Journal  of  Hydrometeorology. 
13:  1631-1638. 

Vivoni,  E.R.  2012.  Spatial  Patterns,  Processes  and  Predictions  in  Ecohydrology: 
Integrating  Technologies  to  Meet  the  Challenge.  Ecohydrology.  5(3):  235-241. 

Vivoni,  E.R.,  Rango,  A.,  Anderson,  C.A.,  Pierini,  N.A.,  Schreiner-McGraw,  A.,  Saripalli, 
S.,  and  Laliberte,  A.S.  2014.  Ecohydrology  with  Unmanned  Aerial  Vehicles.  Ecosphere.  5(10): 
artl30.  http://dx.doi.org/10. 1 890/ES 14-002 17. 1 . 
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Templeton,  R.C.,  Vivoni,  E.R.,  Mendez-Barroso,  L.A.,  Pierini,  N.A.,  Anderson,  C.A., 
Rango,  A.,  Laliberte,  A.S.,  and  Scott,  R.L.  2014.  High-Resolution  Characterization  of  a  Semiarid 
Watershed:  Implications  on  Evapotranspiration  Estimates.  Journal  of  Hydrology.  509:  306-319. 

Pierini,  N.A.,  Vivoni,  E.R.,  Robles-Morua,  A.,  Scott,  R.L.,  and  Nearing,  M.A.  2014. 
Using  Observations  and  a  Distributed  Hydrologic  Model  to  Explore  Runoff  Threshold  Processes 
Linked  with  Mesquite  Encroachment  in  the  Sonoran  Desert.  Water  Resources  Research.  50(10): 
8191-8215. 

Mendez-Barroso,  L.A.,  Vivoni,  E.R.,  Robles-Morua,  A.,  Mascaro,  G.,  Yepez,  E.A., 
Rodriguez,  J.C.,  Watts,  C.J.,  Garatuza-Payan,  J.,  and  Saiz-Hemandez,  J.  2014.  A  Modeling 
Approach  Reveals  Differences  in  Evapotranspiration  and  Its  Partitioning  in  Two  Semiarid 
Ecosystems  in  Northwest  Mexico.  Water  Resources  Research.  50(4):  3229-3252. 

Mascaro,  G.,  Vivoni,  E.R.,  and  Mendez-Barroso,  L.A.  2015.  Hyperresolution  Hydrologic 
Modeling  in  a  Regional  Watershed  and  its  Interpretation  using  Empirical  Orthogonal  Functions. 
Advances  in  Water  Resources.  83:  190-206. 

Browning,  D.M.,  Rango,  A.,  Karl,  J.W.,  Laney,  C.M.,  Vivoni,  E.R.,  and  Tweedie,  C.E. 
2015.  Emerging  Technological  and  Cultural  Shifts  Advancing  Dryland  Research  and 
Management.  Frontiers  in  Ecology  and  the  Environment.  13(1):  52-60. 

Schreiner-McGraw,  A.P.,  Vivoni,  E.R.,  Mascaro,  G.,  and  Franz,  T.E.  2015.  Closing  the 
Water  Balance  with  Cosmic-ray  Soil  Moisture  Measurements  and  Assessing  Their  Spatial 
Variability  within  Two  Semiarid  Watersheds.  Hydrology  and  Earth  System  Sciences.  (In 
Review). 

Anderson,  C.A.,  and  Vivoni,  E.R.  2015.  Land  Surface  States  within  the  Flux  Footprint 
Impact  Land-Atmosphere  Coupling  in  Two  Semiarid  Ecosystems  of  the  Southwestern  U.S. 

Water  Resources  Research.  (In  Review). 

Our  efforts  in  the  last  reporting  period  (August  2014  -  May  2015)  focused  on: 

1)  Continued  data  collection  from  the  environmental  sensor  network  at  the  Jornada 
Experimental  Range  (Fig.  1)  and  at  the  Santa  Rita  Experimental  Range.  Processing  and  analysis 
of  LiDAR  and  UAV  imagery  and  of  the  environmental  sensor  network  data  for  the  use  in 
distributed  hydrologic  models. 

2)  We  summarized  our  analysis  of  the  Cosmic-Ray  Neutron  Soil  moisture  (CRNS) 
measurements  in  a  paper  to  Hydrology  and  Earth  System  Sciences  (Schreiner-McGraw  et  ah, 
2015,  in  review).  Observations  from  the  environmental  sensor  networks  at  both  sites  were  used 
to  validate  the  CRNS  method,  analyze  the  watershed  hydrologic  conditions,  and  improve 
parameterization  of  hydrologic  models.  The  environmental  sensor  network  allowed  us  to 
calculate  the  soil  moisture  over  the  full  watershed  to  validate  the  CRNS  method,  which  is  an 
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important  step  because  previous  validations  had  only  used  point-scale  sensors  that  do  not 
measure  the  same  spatial  scale  as  the  CRNS  method.  We  found  good  correspondence  between 
the  two  methods  (RMSE  =  0.012  and  0.013  m3/m3  at  SRER  and  JER,  shown  in  Fig.  4),  with 
deviations  due  to  bypassing  of  the  CRNS  measurement  depth  during  large  rainfall  events.  This 
limitation,  however,  was  used  in  conjunction  with  measurements  from  the  eddy  covariance 
towers  to  show  that  drier-than-average  conditions  at  SRER  promoted  plant  water  uptake  from 
deeper  layers,  while  the  wetter-than-average  period  at  JER  resulted  in  leakage  towards  deeper 
soils,  as  depicted  in  Fig.  5. 

3)  We  continued  work  on  partitioning  evapotranspiration  (ET)  into  its  components:  bare 
soil  evaporation  (E)  and  plant  transpiration  (T).  This  work  focuses  on  the  Jornada  Experimental 
Range  (JER),  and  was  presented  at  the  Ecological  Society  of  America  100th  annual  meeting 
(Schreiner-McGraw  et  ah,  2015).  We  have  successfully  applied  a  method  to  partition  ET  through 
the  use  of  surface  soil  temperature  measurements  and  the  results  are  promising  (Fig.  6).  We  find 
that  the  T/ET  value  during  the  pre-monsoon  season  is  approximately  0.3  due  to  plant  inactivity. 
When  the  monsoon  starts  and  plants  green  up,  the  T/ET  ratio  shifts  dramatically  and  roughly 
70%  of  ET  comes  from  transpiration.  These  results  are  consistent  with  other  studies  in  the  region 
showing  that  a  simple  method  for  partitioning  ET  is  promising. 

4)  We  expanded  the  environmental  sensor  network  at  the  Santa  Rita  Experimental  Range 
(Fig.  7)  with  the  installation  of  24  soil  sensors  (8  locations  with  sensors  at  3  depths),  4  rain 
gauges  and  1  internal  flume  in  watershed  7  to  complement  the  existing  environmental  sensor 
network  that  was  established  in  watershed  8  in  May  2011.  The  new  measurements  help  quantify 
the  similarities  and/or  differences  observed  in  soil  moisture,  soil  temperature,  rainfall,  and 
internal  runoff  in  the  paired  watersheds,  which  represent  the  landscape  transition  from  desert 
grasslands  (watershed  7)  to  woody  savannas  (watershed  8).  The  measurements  in  both 
watersheds  are  used  in  distributed  hydrologic  models  that  represent  different  woody  plant 
encroachment  densities  depending  on  land  management  decisions. 

5)  We  compared  measurements  from  two  eddy  covariance  towers  at  the  Santa  Rita 
Experimental  Range  separated  by  about  1.5  km  in  distance  and  50  m  in  elevation.  Through  the 
processing  and  analysis  of  the  flux  measurements  over  woody  savanna  landscapes,  we  seek  to 
discern  similarities  and  differences  due  to  land  cover  variations.  As  an  example,  Fig.  8  illustrates 
the  daily  surface  energy  fluxes  at  both  locations  over  a  4  year  study  period. 

These  research  activities  have  demonstrated  a  novel  approach  for  terrestrial  science 
studies  by  integrating  the  deployment  of  environmental  sensor  networks,  UAV  flights  and  HPC 
distributed  hydrologic  modeling  in  a  range  of  sites  in  southwestern  North  America.  Potential 
benefits  from  this  research  include  the  ability  to  deploy  such  sensor  and  modeling  packages  in 
battlefield  environments  for  improved  forecast  of  on-the-ground  conditions,  and  for  the 
environmental  stewardship  of  military  lands  in  arid  and  semiarid  climates. 


4 


^  ^'Legend 


400 

b  Meters 


■  Outlet 
S  ☆  Tower 
S &.•  Raingauges 


.  a  Flumes 


Transect 


Fig.  1.  Chihuahuan  Desert  study  site  at  the  Jornada  Experimental  Range,  New  Mexico,  with 
watershed  boundary  and  stream  network  derived  from  the  1-m  UAV-based  digital  elevation 
model.  The  location  of  the  environmental  sensor  network  sites  are  shown  along  with  two 
photographs  depicting  the  eddy  covariance  tower  and  one  of  the  channel  flumes. 
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Fig.  2.  Example  of  time  series  of  environmental  sensor  network  observations  at  the  Santa  Rita 
Experimental  Range  site  from  June  2011  to  October  2011.  (a)  Basin-averaged  rainfall  rate  from 
5  rain  gauges,  runoff  rate  from  outlet  flume  and  basin-averaged  soil  moisture  at  5  cm,  1 5  cm  and 
30  cm  depth  from  20  profile  locations,  (b)  Basin-averaged  soil  temperature  at  5  cm,  15  cm  and 
30  cm  depth  from  20  profile  locations,  (c)  Latent  and  (d)  sensible  heat  flux  measured  at  the  eddy 
covariance  tower. 
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Fig.  3.  High-resolution  hydrologic  modeling  in  large  watersheds  using  high-performance 
computing.  (a,b)  Comparison  of  spatial  patterns  of  simulated  and  retrieved  (2D-STAR)  surface 
(top  5cm)  Soil  Moisture  (SM)  on  8  and  25  August  2004. (c,d)  Comparison  of  spatial  patterns  of 
simulated  and  estimated  (MODIS)  LST  on  the  same  days.  Model  outputs  are  for  the  overpass 
times  and  are  resampled  at  same  resolution  (800  m  for  2D-STAR  and  1000  m  for  MODIS). 
White  areas  represent  missing  data.  Circles  of  different  size  depict  the  surface  SM  and  LST 
measured  by  individual  stations  at  coincident  times. 
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Fig.  4.  Scatterplots  of  the  spatially-averaged  change  in  soil  moisture  (m3/m3)  derived  from 
CRNS  method  (AGcrs)  and  the  environmental  sensor  network  (A0sn)  at  (a)  SRER  and  (b)  JER. 
The  SEE  and  the  number  of  event  samples  (N)  are  shown  for  each  site. 
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Fig.  5.  Comparison  of  cumulative  fcRNs  and  measured  water  balance  fluxes  (P  and  ET).  CRNS 
estimates  of  infiltration  (I),  outflow  (O)  and  leakage  (L)  are  either  depicted  as  cumulative  fluxes 
(O  =  ET  +  L)  or  as  total  amounts  during  the  study  period  (I  and  L)  as  arrows  in  the  soil  water 
balance  box  of  depth  z*.  Shaded  regions  indicate  the  summer  seasons  (July-September).  The 
horizontal  line  represents  fcRNs  =  0. 
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Fig.  6.  Evapotranspiration  partitioning  at  the  JER  for  the  summer  of  2013. 
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Fig.  7.  Location  of  the  study  watersheds  7  and  8  (WS  7  and  WS  8)  within  the  Santa  Rita 
Experimental  Range  (SRER).  The  0.31  meter  aerial  photographs  are  products  of  the  LiDAR 
flight  taken  in  April  2011.  The  environmental  sensor  network  was  expanded  into  WS  7  with  the 
addition  of  24  soil  sensors  (at  3  depths  at  8  locations)  to  measure  soil  moisture  and  temperature, 
4  rain  gauges  to  measure  precipitation,  and  1  internal  mini-flume  to  measure  runoff. 


8 


Rainfall  [mm] 


Daily  Average 


Jan12  Jan13  Jan14  Jan15 


Jan12 


Jan13 


Jan14 


Jan15 


Fig.  8.  Processed  fluxes  measured  at  the  two  eddy  covariance  towers  (ASU  and  collaborating 
USDA-ARS  measurements),  including  (a)  sensible  heat  flux,  (b)  latent  heat  flux,  and  (c)  ground 
heat  flux,  at  a  daily  averaged  time  scale. 
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